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Abstract 

Myers-Perry black holes with a single spin in d > 5 have been shown to be unstable if 
rotating sufficiently rapidly. We extend the numerical analysis which allowed for that result 
to the asymptotically AdS case. We determine numerically the stationary perturbations that 
mark the onset of the instabilities for the modes that preserve the rotational symmetries of the 
background. The parameter space of solutions is thoroughly analysed, and the onset of the 
instabilities is obtained as a function of the cosmological constant. Each of these perturbations 
has been conjectured to represent a bifurcation point to a new phase of stationary AdS black 
holes, and this is consistent with our results. 



1 Introduction 



There has been recent progress in understanding the phase diagram of higher-dimensional asymp- 
totically flat vacuum black holes. The purpose of this work is to extend some of the techniques 
used to the case of asymptotically anti-de Sitter (AdS) black holes. 

Let us first review the asymptotically flat case. While the four-dimensional Kerr black hole is 
unique and (expected to be) stable, the higher-dimensional picture is much richer pQ. In d = 5, the 
discovery of the black ring by Emparan and Reall [2j showed that it can have the same conserved 
charges as the Myers- Perry (MP) black hole [3J. This non- uniqueness discovery triggered the 
research in higher dimensional black holes and a variety of explicit rotating solutions have been 
recently found: black Saturns [1], concentric rings 016], orthogonal rings 018], and generalisations 
thereof. Consider solutions rotating on a single plane which are in thermal equilibrium (i.e. in the 
case that there are disconnected components of the event horizon, these have the same temperature 
and angular velocity). One interesting fact is that in the zero temperature limit the new solutions 
and the MP black hole coincide in the same nakedly singular solution [9j. In d > 5, this limit does 
not exist, and the MP black hole has an unbounded angular momentum for a given mass. In a 
certain sense, the singular limit is 'resolved'. However, does the MP black hole still connect to the 
new solutions? In d > 5, the black ring and the other solutions have not been constructed exactly, 
but are expected to exist. In fact, approximate methods for solving the Einstein equations in the 
ultraspinning limit (large angular momentum with respect to the mass scale) indicate that they do 

exist (ini nn us] . 

The answer seems to be that the MP black hole connects to the other solutions through yet 
more families of black holes. These families bifurcate from the MP branch at the (stationary) onset 
of the so called 'ultraspinning instabilities'. The first step in understanding this connection was 
given by Emparan and Myers [13], who argued that the singly-spinning (i.e. rotating on a single 
plane) MP black hole should be unstable if rotating sufficiently rapidly. The reasoning is that if 
the angular momentum per unit mass is large enough, MP black holes start behaving like black 
branes, since their horizons become disk-like along the rotation plane in this ultraspinning limit. 
But it has been shown by Gregory and Lanamme [TS] that black branes are unstable and therefore 
one concludes that MP black holes should be unstable for a sufficiently large (but finite) angular 
momentum per unit mass. Emparan and Myers further pointed out that at the onset of such an 
instability, there should exist a stationary perturbation that preserves the rotational symmetries 
of the background. This perturbation would then signal the existence of a new family of black 
holes with spherical horizon topology. This is analogous to what happens in the black brane case: 
the threshold mode of the Gregory-Laflamme instability is the perturbative signal for a family of 
non- uniform branes |X6(, 117]. Ref. [10] conjectured that the new families bifurcating from the MP 
branch, when continued along the phase diagram, would connect continuously to the black ring 
branch, the black Saturn branch, and so on. Therefore, the connection of non-uniqueness with 
instabilities provides a partial understanding of the proliferation of higher dimensional black hole 
phases. 

Recently, these conjectures were put on a firmer footing when it was shown numerically that 
the singly-spinning MP black hole does indeed possess such threshold modes pS [H] . Further- 
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more, these perturbations exhibit an underlying harmonic structure and induce deformations on 
the shape of the horizon which are consistent with the proposal that the new families will connect 
to the black ring, to the black Saturn, etc. Ref. |18] also conjectured that these classical insta- 
bilities can only occur when the black hole possesses at least two distinct local thermodynamic 
instabilities, as we shall review later. This is a necessary but not sufficient condition. The reason 
is that local thermodynamic instabilities are associated to the lowest harmonics, and thus with 
the asymptotic charges, the mass and the angular momentum. Higher harmonics, which cannot 
change the asymptotic charges and are associated to the bifurcation to new families, should only 
become unstable for rotations higher than the thermodynamic modes. This conjecture allows for 
any number of independent angular momenta, and has already been verified in the instability of 
d > 5 cohomogeneity-1 MP solutions (equal independent spins in odd d) |20j . This sector of MP 
black holes has a regular extremal limit and the angular momentum is bounded from above, yet the 
conjecture of Ref. [18] indicates that an instability is possible. Indeed linear perturbations growing 
exponentially with time were found, which is an important check since the previous work could 
only determine the stationary threshold modes. 

In this paper, we will extend the results in Refs. [181 IS] to the asymptotically AdS case. 
The properties of asymptotically AdS spacetimes have been greatly explored due to the AdS / CFT 
correspondence, which equates quantum gravity (string theory) in an asymptotically AdS spacetime 
to a conformal field theory (CFT) living on the boundary of that spacetime |21[ I22|. In particular, 
the phase diagram of black holes in AdS is in direct correspondence to the phases of the dual CFT 
at finite temperature. Therefore, even if many higher-dimensional black hole solutions present 
classical instabilities, or do not dominate thermodinamically the gravitational partition function, 
they still provide a valuable insight into the CFT phases. 



Black rings and black Saturns in AdS have been constructed in certain approximations [23J2I and 

the conjectures of Ref. [10] regarding the connections between singly-spinning black hole families 
were extended to the asymptotically AdS case also in Ref. [23]. Notice that singly-spinning MP- 
AdS black holes [25J (and indeed any asymptotically AdS stationary black hole |26j ) present a 
BPS-type upper bound on their angular momentum, \J\ < Mi (with M being the mass and t the 
cosmological length), while there is no such bound in the asymptotically flat case. However, the 
instabilities/bifurcations are still expected as we increase the rotation up to that bound [23] (see 
Fig. In this paper we will numerically determine the dependence of the stationary threshold 
modes on the cosmological constant. The procedure is an extension to d > 6 of the d = 4 Kerr- 
AdS analysis [2Z], where no classical instability was found; the latter was itself an extension of the 
Schwarzschild-AdS case |28j. 

We confirm the conjecture of Ref. |18j . which gives a necessary but not sufficient condition for 
the onset of the ultraspinning instability. In this paper we generalise the previous formulation so 
that it can be applied to the asymptotically AdS case. We find that all bifurcations occur when 
the MP- AdS black hole is unstable under superradiance [2H1 1^1 EI] , which implies that the new 
families of black holes will inherit this instability. 

Before proceeding, let us mention that asymptotically flat singly-spinning MP black holes also 

1 See Ref. [21] for the extension of the more systematic blackfold approach of Refs. [Ill 1121 IT5] to asymptotically 
AdS black holes. 
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Figure 1: Phase diagram of singly-spinning MP-AdS black holes in d > 6. We plot the entropy 
S vs. the angular momentum J, at a fixed value of the mass M, in units of the AdS curvature 
radius t. The figure illustrates the conjecture of Ref. |23j . At sufficiently large spin the MP- AdS 
solution becomes unstable for axisymmetric perturbations (dashed line), and at the threshold of 
the instability a new branch of black holes with a central pinch appear (/). As the spin grows, new 
branches of black holes with further axisymmetric pinches (II, III, ■ ■ ■ ) appear. We determine 
numerically the points where the new branches appear, but it is not yet known in which directions 
they run. 

suffer from an instability that breaks their axisymmetry [32, 33 1. This occurs even in d = 5, and for 
d > 6 it sets in for lower rotations than the ultraspinning instability that we study here. However, 
the threshold mode is not associated to new stationary black holes. It would be interesting to find 
whether this non-axisymmetric instability also extends to asymptotically AdS black holes and to 
understand how it relates to the AdS superradiant instability. 

This paper is organised as follows. In section [2j we review the properties of singly-spinning 
MP- AdS black holes and discuss the ultraspinning instability. In section [3j we detail the linear 
perturbations problem that is solved numerically. We conclude in section |4j with the discussion of 
the results. 

2 Myers-Perry- AdS black holes 

In this section, we will start by reviewing the properties of singly-spinning MP-AdS black holes. 
Although for such solution the angular momentum (at fixed mass) has an upper bound, we provide 
a thermodynamic argument, discussed in subsection |2.2[ according to which these black holes may 
be afflicted by the ultraspinning instability. This motivates the numerical work, described in later 
sections, of searching for the onset of the instability. We conclude this section by noting that all 
the AdS black holes potentially afflicted from the ultraspinning instability are also afflicted by the 
AdS superradiant instability. 
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2.1 Solution 



The four-dimensional Kerr-AdS black hole was found by Carter |34j . and its extension to higher 
dimensions, when rotation on a single plane is considered, was obtained by Hawking, Hunter and 
Taylor |25j The metric is given by 



ds 2 
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where dQ 2 d _^ is the line element of a unit-radius (d — 4)-sphere and 
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A e = r 2 + a 2 cos 2 9 , £ = r 2 + a 2 cos 2 9 , E = l-^. 

(2.1) 

This solution of the Einstein equations with negative cosmological constant, = —(d— 1)£~ 2 g^ u , 
is parameterised by three length scales, namely the AdS curvature radius I, the mass-radius r m 
and the rotation parameter a. The Komar mass M and angular momentum J of the black hole are 
given in terms of these parameters by 



AT 



A 



d-2r[ 



d-3 



8-kGE 2 



1 + 



d-4 



J 



A, 



d-2 r m 3 „ 



8ttG H 2 



(2.2) 



where Ad-2 = 2 7r^ 1 ^ //2 /r[(d — l)/2] is the volume of a unit-radius (d — 2)-sphere, and G denotes 
Newton's constant which we set to one (G = 1). The solution satisfies the BPS-like bound [26J 



\J\ < Mi 



a < 



(2.3) 



The metric (2.1) does not describe a black hole if this bound is saturated. In the limit \a\ —> £ 
either the charges (mass and angular momentum) and the horizon equator circumference of the 
black hole diverge, or the horizon vanishes and the solution is nakedly singular. 

The event horizon lies at the largest real root r = r + of A r (r) = 0. The horizon angular velocity 
measured with respect to a non-rotating frame at infinity is [351 136] 
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(2.4) 



while the temperature and entropy of the black hole are [35] 



S = 



1 



4tt r_i-(a 2 



d- 3 + (d- 1) 



„2 1 
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(2.5) 



Since the Hawking-Hunter- Taylor solution usually refers to the five-dimensional black hole with two independent 
spins, also found in [23], we will refer to the black hole represented by (2.1 1 simply as (singly-spinning) Myers-Perry- 
AdS black hole. 
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In d = 4, the temperature vanishes for \a\ = r +\J (^ r + + £ 2 )(£ 2 — r^) -1 . However, as a con- 
sequence of the constraint ( |2.3[ ), we have an extremal regular (i.e. with finite size horizon area) 
black hole if and only if r + /£ < In d = 5, like in the asymptotically flat space limit £ — > oo, 

the temperature goes to zero at the singular limit \a\ —> r m , whereas the bound on the rotation 



parameter for a finite r + is (2.3). For d > 6, the temperature can never vanish and thus there is no 



extremal solution. The parameter space is bounded only by (2.3). This is in sharp contrast with 
the asymptotically flat case, in which MP black holes rotating on a single plane in d > 6 have no 
upper bound on their angular momentum. 

2.2 Thermodynamic zero-modes and the ultraspinning regime 

The condition of local thermodynamic stability for a black hole with charges M, J{ and entropy S 
is the positivity of the Hessian 

We write the condition in its most general form, i.e. taking into account the several angular momenta 
allowed in higher dimensions, i = 1, . . . , \_(d — l)/2j (where |_-J stands for the smallest integer 
part). In this section, we will review the arguments of Refs. [18, 20J which relate the eigenvalues 
of this thermodynamic Hessian to the occurrence of classical instabilities, namely the ultraspinning 
instability. 

There are two general properties of this Hessian for AdS black holes: (i) For all (non-extremal) 
asymptotically flat vacuum black holes, in any d, the Hessian has at least one negative eigenvalue 
|20j ; so, by continuity, black holes whose size is much smaller than the AdS curvature radius 
should be thermodynamically unstable, (ii) Large black holes in AdS, however, are expected to be 
thermodynamically stable, and thus the Hessian should be positive definite. 

Let us take the d = 4 Kerr-AdS case, where we have a 2 x 2 thermodynamic Hessian. One of 
the eigenvalues is positive in the entire parameter space. However, the second eigenvalue can be 
negative, null or positive depending on the parameters. The null point is where the specific heat at 
constant angular velocity diverges, changing sign. We shall say that, at the surface in parameter 
space where the Hessian becomes degenerate, we have a thermodynamic zero-mode, corresponding 
to the eigenvector with zero eigenvalue. The situation is analogous in the singly-spinning d = 5 
case (say J2 = 0), with two eigenvalues being always positive, and one changing sign. In the d > 6 
singly-spinning case (say Jj = for i > 1) things are more interesting. Two of the eigenvalues of the 
thermodynamic Hessian change sign, and all the others remain positive throughout the parameter 
space. So there will be two thermodynamic zero-modes, marking the onsets of two distinct local 
thermodynamic instabilities. 

Let us consider the parameterisation (r+ , a) , which uniquely specifies a singly-spinning MP- AdS 
black hole, for a given AdS curvature radius I. Recall that \a\ < £ is the bound on the parameter 
space, as represented in Figure [2] for d = 6 (it is qualitatively similar in d > 6). Now, for large r+/£, 
the black hole is thermodynamically stable. As we decrease r+, one of the eigenvalues changes sign 
first, which corresponds to the transition between regions A and B in Figure [2j Let us label the 
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Figure 2: Parameter space of singly-spinning Myers-Perry— AdS black holes in d = 6: dimensionless 
rotation parameter a/r m as a function of the dimensionless mass-radius parameter r m jt. Regular 
black holes exist only for a < i. The thermodynamic I = and I = 1 zero-mode curves, described 
by (2.7) and (2.8), are plotted. The thermodynamic I = 1 zero- mode curve defines also the 
ultraspinning surface, above which (region C) the black holes might be ultraspinning unstable 
(later section [4] and Fig. [4] confirm this is indeed the case) . The superradiant curve = 1 is 
also plotted. Above this curve all black holes are superradiant unstable. For d > 6 the plot is 
qualitatively very similar. 



associated zero-mode as the I = zero-mode, for reasons that will become apparent later. Its 
parameter space locus isj^] 

d 3 

Thermodynamic I = zero-mode (d > 4): r\ = — — — (a 2 + I 2 + yV - -f d a 2 £ 2 + ¥) , (2.7) 

where ^ = 2 (d 2 - 6d + l) /(d - 3) 2 . 

The I = zero-mode is present for any d > 4. However, for d > 6 there is an additional zero- 
mode, corresponding to the second eigenvalue that becomes negative as we lower r + for fixed a. 
This corresponds to the transition between regions B and C in Figure [2} Let us label the associated 
zero- mode with 1 = 1. Its parameter space locus is: 

Thermodynamic I = 1 zero-mode (d > 6): r\ = —- - *j- (a 2 + I 2 - - Id a 2 £ 2 + . (2.8) 

What happens as we decrease r + further for fixed a? As discussed in the Introduction, singly- 
spinning MP black holes were conjectured to be unstable if rotating too rapidly, i.e. \a\ S> r + , 



Notice that the Hawking-Page transition [37] occurs for r+/£ larger than (2.7 1. This implies that hot AdS space 
is preferred (in the grand-canonical ensemble) over these black holes if r+ < £, in any d. Global thermodynamic 
stability is stricter than local thermodynamic stability. 
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since they start behaving like black branes [14] . In the asymptotically AdS case, there should also 
be an analogous unstable regime [23J. Now, there can be two types of instabilities. One for which 
the perturbations break the rotational symmetry of the background, i.e. the ds Killing vector of 



(2.1); see \32\ I33j for the analysis of these perturbations in the asymptotically flat case. And one 
for which the spatial isometries of the background are preserved. Instabilities of the latter type 
admit stationary threshold modes, which may correspond to bifurcations to new black hole phases 
[T4"l [TO] . Ref . [2] suggested that an order of magnitude estimate for the critical rotation at the 
bifurcation point could be given by what we described above as the I = 1 zero mode. In fact, this 
point marks the transition to the black brane-like behaviour and it only occurs for d > 6. The 
subsequent understanding of the problem allowed for a more precise formulation of the conjecture, 
which can also be extended beyond the singly-spinning case. 

Indeed, the ultraspinning conjecture of Ref. [18] proposes that instabilities whose onset is a 
stationary and axisymmetric mode, i.e. a mode which gives a bifurcation to a new stationary black 
hole family, can only occur after the I = and I = 1 zero-modes, i.e. in region C of Figure [2] in the 
present case. This is a necessary but not sufficient condition. The reasoning is that the onset of 
these instabilities corresponds to the I > 2 zero-modes, as we shall review below. 

Ref. }38J, in the static case, and Refs. [391 120]) in the general case, showed that a zero- mode 



(or a negative mode, the eigenvector of a negative eigenvalue) of the thermodynamic Hessian (2.6) 
is also a stationary and axisymmetric zero-mode (negative mode) of the black hole Euclidean 
action. Since we are dealing with stationary and axisymmetric modes, this result extends to 
the Lorentzian action. In the case of a zero-mode, it corresponds to a classical perturbation of 
the black hole which preserves its temperature and angular velocities, but changes its asymptotic 
charges (mass and angular momenta) . The change in the asymptotic charges can simply be inferred 



from the eigenvector (5M, 5Ji) of the Hessian (2.6) associated with the null eigenvalue. In the case 



of a negative mode, it corresponds to both: (i) an off-shell perturbation of the black hole which 
makes the (Euclidean) gravitational partition function pathological, the expected signal of a local 
thermodynamic instability; and (ii) an on-shell stationary perturbation of a black brane, where 
the negative eigenvalue corresponds to (minus) the 'mass-squared' coming from the dimensional 
reduction along the brane directions. This connection is the basis of the Gubser-Mitra conjecture 



The first example of the relation between negative modes of the action and thermodynamic 
stability was the determination of the negative mode of the Schwarzschild black hole [UJ. The 
existence of such mode is a direct consequence of the fact the Schwarzschild black hole has negative 
specific heat. Indeed, the extension to the Schwarzschild- AdS case [28] showed that the relevant 
mode changes sign exactly where the specific heat changes sign, as expected. This result was 
extended to the Kerr- AdS case [27] , where the numerical techniques used in [181 ES EH 02] , and 
also in this paper, were applied first to a black hole stability problem. 

Now we can justify the labels I = and I = 1 used before, following [18J. They are based on 
the harmonic structure of the metric perturbations corresponding to the zero- modes. The I = 
zero-mode has no nodes, i.e. the metric perturbation does not vanish anywhere on the horizon, 
whereas the I = 1 zero-mode has one node, as plotted in [18] for the singly-spinning MP case. Now, 
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the I = and I = 1 zero-modes are associated with the asymptotic charges of the spacetimej^] This 
is consistent with the fact that they are determined from the thermodynamic Hessian ( |2.6[ ), the 
changes in the charges being proportional to the eigenvector (SM, <5Jj). Notice that, being signalled 
by the Hessian of S(M, Jj), i.e. by the equation of state of the background black hole family, the 
I = and / = 1 zero-modes can only change the black hole into another black hole of the same 
family, e.g. another MP black holej^] 

What is found by solving the problem of axisymmetric perturbations is that instabilities appear 
with new stationary zero-modes, not predicted by the thermodynamic Hessian. We refer to them 
as non-thermodynamic zero-modes. Their harmonic structure - the number of nodes, found to 
be I > 2 - is consistent with the deformations of the event horizon proposed in Refs. |1Q|. 123]: see 
Figure[l} These higher harmonics cannot change the asymptotic charges, as exemplified explicitly in 
|20j . The perturbed black holes will have the same mass and angular momenta as the unperturbed 
solutions, and also, as imposed by the boundary conditions, the same temperature and angular 
velocities. Therefore, either they correspond to trivial gauge modes or to bifurcations to new black 
hole solutions. We have shown in specific cases |20|. 119] that no regular gauge modes are allowed. We 
conclude that the I > 2 zero-modes mark a black hole bifurcation, and also the onset of a classical 
ultraspinning instability, one per each mode. That these linear modes do grow exponentially with 
time was verified in |20j. 

The ultraspinning conjecture of Ref. [TH] predicts that zero- modes with I > 2 can only occur 
for rotations higher than that of the I = 1 zero- mode. In the case of singly-spinning MP-AdS black 
holes, this corresponds to the region C of Figure [2] This is based on the examples known, and on 
the intuition that higher harmonics should be more stable than lower harmonics. In general, when 
several distinct angular momenta are present, several different zero- modes with / = 1 may exist, 
as exemplified in |20j . So the ultraspinning surface (in our case, the surface between B and C) 
is defined as the surface in parameter space which encloses the region where the thermodynamic 
Hessian (2.6) has less than two negative eigenvalues (in our case, regions A and B). 

Since this thermodynamic criterion is a necessary but not sufficient condition, to confirm the 
presence of these instabilities we still have to solve the linear perturbations problem, and look for 
the I > 2 non-thermodynamic zero-modes. This study is done in section [3} 



2.3 Superradiant instability 

Before studying the ultraspinning instability of the singly-spinning MP-AdS black hole, we recall 
that this is not the only instability present in the system. Indeed, it competes with the well-known 

4 Ref. [20] analysed MP black holes with equal spins in odd d, which have the convenient property of being 
cohomogeneity-1. In that case, the number I labels a precise harmonic of the CP N base space, and it can be shown 
that the I — harmonic is associated to the mass only (indeed <5J; = in the Hessian eigenvector), while the 
I — 1 harmonic is associated to the angular momenta only (indeed SM = in the Hessian eigenvector). In the 
singly-spinning black hole both the I = and I = 1 'harmonics' change the mass and the angular momentum. 

5 In the present paper, these zero-modes give another singly-spinning MP-AdS black hole. However, in the 
cohomogeneity-1 MP(-AdS) case [20], the I = 1 zero-mode breaks the symmetries between the spins, as predicted by 
the eigenvector (5M,SJi) having distinct SJi. The zero-mode therefore takes the black hole out of the equal spins 
sector, but still in the general MP(-AdS) family. 
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superradiant instability that afflicts rotating AdS black holes whenever Slut- > 1 [Z2 [30], 131] . 
The superradiant instability is associated to perturbations that break the symmetry generated 
by d<j), i.e. it requires perturbations whose angular momentum quantum number m along the <j>- 
direction is non-vanishing. The mechanism of the instability is simple. A gravitational wave 
co-rotating with the background black hole can extract rotational energy from the black hole if its 
frequency u satisfies the relation uj < raOg, where Cljj is the angular velocity of the black hole. 
This phenomenon is known as superradiance. On the other hand, the AdS gravitational potential 
effectively acts as an asymptotically reflecting box. The multiple superradiant amplifications and 
reflections drive the system unstable if and only if Q,h& > 1- 

Since there is a competition between the two instabilities, it is relevant to ask whether the 
ultraspinning instability occurs when Qjji < 1. We find a negative answer. Indeed, in Fig. [2] we 
represent the superradiant surface, £Ih@ = 1, by a dashed line. It lays entirely in region B. Black 
holes above this surface have f2#^ > 1, which includes all the black holes in region C that might 
also be afflicted by the ultraspinning instability. Therefore, all the AdS black holes that suffer from 
the ultraspinning instability that we identify in the next section are also likely to be unstable under 
superradiance. 



3 Perturbations of MP- AdS black holes 

In this section, we will detail the linear perturbation problem, which consists in studying the 
spectrum of the Lichnerowicz operator on the black hole background. We will list the equations to 
be solved and the boundary conditions, and we will briefly discuss their numerical implementation. 



3.1 The Lichnerowicz eigenvalue problem 

To find the onset of the ultraspinning instability we have to solve the linearised Einstein equations 
with a negative cosmological constant. We search for stationary modes that depend on the radial 
and polar coordinates, r and 9, but preserve the R x U(l) x SO(d — 3) symmetries of the back- 
ground MP- AdS black hole (2.1). The general ansatz for the perturbed metric h^ u satisfying these 
conditions is|3 



+ 



A r (r) 
E(r,0) 

em; 

A r (r) 



dt 



a sin 



■ sin 2 9 Aq(0) Sui 
E(r,0) 



r 2 + a 2 



ae~ Suj dt 



e s ^° [dr - 5 X sin 9 d9} 2 + §^ e^ 1 d9 2 + r 2 cos 2 9 e 5 * dQ 2 d _ 4) 

Ag{9) v ; 



(3.1) 



where {Suq, Sv%, Sfio, 6u, 6x, 8$} are small quantities that describe our perturbations, and are 
functions of (r, 9) only. We will solve numerically the coupled partial differential equations (PDEs) 



Stationarity and axisymmetry lead to a metric perturbation independent of t and (f>, respectively. The transverse 
S d ~ 4 line element is left unperturbed, apart from its (r, #)-dependent scale factor, because we wish to preserve the 
SO(d — 3) isometry. On the other hand, the components tr, t9, cjyr, <f>9 can be discarded: the corresponding TT gauge 
conditions and the linearised equations of motion can be solved by separation of variables in r and 9, and regularity 
implies that those components vanish. 
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that govern these perturbations. We choose to work in the traceless-transverse (TT) gauge, 



^ = and V M V = • ( 3 - 2 ) 
The linearised Einstein-AdS equations are: 

(A L h)^ = -V p V p h^ - 2 Rfjhpa = , (3.3) 
where the operator A^, relates to the standard Lichnerowicz operator Al as 

(A L h)^ = (A L h)^ + 2{d - l)r 2 V • ( 3 - 4 ) 

Following the strategy used to study the ultraspinning instability in asymptotically flat black holes 
in [2"0"1 ITU] , we will actually consider a more general eigenvalue problem, 

(A L h)^ = Xh^. (3.5) 

More concretely, we will be looking for negative modes (A < 0) of the operator A^. 

This eigenvalue problem arises in two instances. One is the computation of quadratic quantum 
corrections to the gravitational partition function in the saddle point approximation [UJ (see |27j for 
the application to the Kerr-AdS black hole). The quantum corrections present a pathology whenever 
there is a negative mode. As we have mentioned, each negative eigenvalue of the thermodynamic 



Hessian (2.6) gives a negative mode, which we refer to as being thermodynamic. 



In the asymptotically flat case, (3.5) also governs the gravitational perturbations of the form 



of the rotating black string that is constructed by adding a flat extra dimension z to 



the MP geometry. The AdS analogue are the warped AdS black strings constructed in Ref. 
These solutions are the generalisations of the standard warped AdS black string, whose stability 
was studied in |44j . in the case where the transverse black hole is also asymptotically AdSQ 

Our strategy (following Refs. [271 H]) to study the classical perturbations (A = 0) of the 



black hole is to look for a solution of (3.5), i.e. a negative mode of the black hole, and then 
vary the rotation parameter a of the black hole until the negative mode becomes a zero-mode, 
i.e. until A — > 0. This strategy is motivated by the availability of powerful numerical methods 



for solving eigenvalue equations of the form (3.5). This strategy is motivated by the availability 



of powerful numerical methods for solving eigenvalue equations of the form (3.5). In particular, 
after an adequate discretisation scheme using spectral methods [16], (3.5) reduces to an algebraic 
eigenvalue problem which we then solve using the inbuilt routine Eingenvalues of Mathematica. 

3.2 Boundary conditions 



To solve the Lichnerowicz eigenvalue system of equations (3.5) we have to impose boundary condi- 



tions on the metric perturbations (3.1). More concretely, we have to specify boundary conditions 
at the horizon, r = r + , at the asymptotic infinity, r — > oo, at the rotation axis 9 = 0, and at the 
equator 8 = ir/2. The strategy to find the appropriate boundary conditions was already discussed 



7 These should not be confused with the uniform black string solutions of Ref. 45 , whose stability problem is not 
given by Eq. (|3.5[). 
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in great detail in section 4 of Ref . |19j , where the ultraspinning instability in asymptotically flat MP 
black holes was studied. The discussion there translates straightforwardly to the present situation, 
and hence we shall be brief. 

We impose regularity of the metric perturbations on the event horizon by demanding regu- 
larity of the Euclideanised perturbed geometry. The idea is to peform the standard Euclidean 
continuation of the black hole metric, and find the conditions required to have regularity at 
the bolt (Euclideanised horizon). For the background black hole, this demands that we iden- 
tify (r, 0) ~ (t, (f) + 2tt) ~ (r + /3, — if2#/3), where r = it is the Euclidean time, (3 = 1/Th is the 



inverse of the black hole temperature (2.5), and Q,h is its angular velocity (2.4). The boundary 



conditions for the metric perturbations h^ v can now be determined by demanding that h fiu dx^dx l/ 
is a regular symmetric 2-tensor on the background manifold. This requires the following boundary 
conditions at the horizon, located at p = [4(r — r + )/A^,(r + )] 1 / 2 = 0, 



5x, Suj = 0(p 2 ) , 5uq - 5p = 0(p 2 ) 



d p 8pi, d p 5vi, d p 5<5> = O(p) . 



(3.6) 



Note that we have not just imposed regularity of the perturbed metric. If hp V dx^dx v is a regular 



2-tensor, then perturbations obeying (3.6) preserve the angular velocity and temperature of the 



background black hole. This is a fundamental property to make the connection between the classical 
ultraspinning instability and black hole thermodynamics that was discussed in section [2] [20J. 

To find the boundary conditions that the metric perturbations must satisfy at the axis of 
rotation (0 = 0), where vanishes, we proceed in the same way and require that h^dx^dx" is 
a regular symmetric 2-tensor, i.e. the components hp U are regular when expressed in coordinates 
where the background metric components are regular. Regularity as — > demands: 



5ui - 5pi = O(0 2 ) , dgSx, d e 5u, d e 5p , <%<5f , d e 5® = 0(6) . 
Analogously, at the equator, 9 = ir/2, we impose the boundary conditions: 

5x = O(x) , 5$ - = 0(x 2 ) , d x 5u, d x 5p , d x 5v 0) d x bv x = 0(x) 
as x = cos 8 — > 0. 



(3.7) 



(3.8) 



We have explicitly checked that the boundary conditions (3.6), (3.7) and (3.8) are consistent 



both with the Lichnerowicz eigenvalue equations (3.5) and the TT gauge conditions (3.2). Indeed, 



the first term in the series expansion of the eigenvalue equations vanishes after we impose (3.6), (3.7) 



and (3.8). On the other hand, we can use the TT gauge conditions to express, e.g. {5vq, 8v\, 5$} as 
functions of {5po, 5p±, 5lj, 5x} and their first derivatives. Again, the first term of a series expansion 



of these TT gauge conditions is consistent with (3.6), (3.7) and (3.8). 

At spatial infinity, r — >• oo, our second order equations of motion allow for two radial depen- 
dences of the metric perturbation. We require that the boundary conditions preserve the asymp- 
totics of the background spacetime, in the sense that not only (in our Boyer-Lindquist coordinates) 

1 



h 



^ u \r— s>oo 



o, 



(3.9) 



for some constant a > that depends on the particular metric component and the number of 
spacetime dimensions d, but also in the sense that the perturbations are normalisable. The negative 
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modes are "gravitons with a positive mass" , and normalisability requires that we pick the fastest 
decaying mode (cf. the non-rotating case [28] )■ This happens to be the mode that our numerical 
code can identify straightforwardly. 



3.3 Imposing the TT gauge conditions and the boundary conditions 



We wish to solve the Lichnerowicz eigenvalue problem (3.5) for the seven metric perturbations 



described in (3.1), namely {6/j,q, 8jj,\, 5\, 8oj, 8uq, 8u\, 5&}, subject to the TT gauge conditions (3.2). 



The latter eliminate three functions. We choose to solve the gauge conditions (3.2 ) for {<5fo, 8v\, 5$} 



in terms of {5/j.o, Sfii, 5x, 8lo} and their first derivatives. Plugging this information in the full set 
of the perturbation equations ( |3.5[ ), we find that only four equations remain of second order in 
{dfio, 6(j,i, 6x, 8u}. Explicitly, these equations are: 

(Ax/t) rr = A h rr , 
{A L h) r6 = \ h r9 , 
(A L h) ee = Xhee, 



a(A L h) tt + 



(r 2 + a 2 + a 2 sin 2 9) 



(r 2 + a 2 ) sin 2 9 



[A L h) t4 , + 



(r 2 + a 2 ) sin 2 9 



(Ai 



(3.10) 



A 



a h tt + 



E (r 2 + a 2 + a 2 sin 2 9) 
(r 2 + a 2 ) sin 2 9 



ht(f> + 



(r 2 + a 2 ) sin 2 

and they describe the final set of equations that we have to solve. A non-trivial consistency check of 



our procedure is to verify that the final equations (3.10) imply that the remaining equations in (3.5) 



are automatically satisfied (the latter equations are of third order in the independent perturbation 
functions once the TT gauge conditions have been imposed). We have explicitly verified that this 
is the case. 



We will solve numerically the final eigenvalue problem (3.10) using spectral methods (see e.g. 
). To do so, we find it convenient to introduce new radial and polar coordinates, 

r 



y 



l 



cos 9 , 



(3.11) 



such that < y < oo and < x < 1 . The implementation of the method is simpler if all functions 
obey Dirichlet boundary conditions on all the boundaries. Therefore, we redefine our independent 
functions according to: 



qi(y,x) 
Q3(y,x) 



r 

r+ 
r 



11 x{\ 
II x(l 



x)5fx (y,x) , 
x)8fii(y,x) , 



q2(y,x) = r m (l -x)Sx(y,x) 
qi{y,x) = x5u(y,x) , 



(3.12) 



so that the q^s vanish at the boundaries. This guarantees that the boundary conditions (3.6), (3.7) 



(3.8) and (3.9) are correctly imposed. 

In the limit £ — > oo, our equations and results reduce to those of the asymptotically flat case 
studied in [THl[Tj5]. In this case, we explicitly proved that the ultraspinning zero-modes that we find 
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cannot be pure gauge modes. That is, there is no pure gauge perturbation, obeying the boundary 
conditions we impose, which could potentially generate the regular metric perturbations that we 
consider. By continuity, we expect this property to hold true when the cosmological constant is 
switched on. 



4 Results and discussion 

The problem of finding the onset of the axisymmetric ultraspinning instability in the singly-spinning 



MP-AdS black hole reduces, at this point, to solving the equations (3.10) for the metric perturba- 



tions (3.12) in the TT gauge and that obey Dirichlet boundary conditions. The outcome of this 
analysis are the dimensionless negative modes, Ar^, of the modified Lichnerowicz operator A^,. 
These are obtained for each pair of parameters that specify the black hole for a fixed AdS radius 
£, namely the dimensionless mass radius, r m /£, and the dimensionless rotation parameter, a/r m . 
We will later describe the results in terms of the dimensionless asymptotic charges M/£ d ~ 5 and 
J/£ d ~ 2 , which are more physical quantities. However, r m and a, which appear in the metric (|2.1l), 



are more convenient for the numerical analysis. 

Examples of the spectrum of negative modes in d = 6, as a function of the dimensionless rotation 
parameter, for two particular values of r m /£, are presented in Fig. [3j We expect that for higher d 
the spectrum is qualitatively similar to d = 6 (we explicitly confirmed this for d = 7). 

The first property to notice is encoded in the curve on the left plot and on its counterpart on 
the right plot, which is the first curve counting from the left. This curve describes what we called 



the thermodynamic I = negative mode in subsection 2.2 As discussed there, it was proven in 
Ref. [20j that this negative mode is always present if the black hole is asymptotically flat. However, 
the AdS curvature radius introduces a new scale in the system that changes this property. For black 
holes with small r m /£ (namely for r m /£ < 0.906 in d = 6) this negative mode is present for any 
value of a/r m , see e.g. T m /£ — 0.5 in the right plot of Fig. [3j while for black holes with large v m /£ 
(namely for r m /£ > 1.587 in d = 6) the thermodynamic / = eigenvalue is positive, independently 
of the rotation. When we fix r m /£ to be in between these two regimes, the thermodynamic / = 
eigenvalue is positive for small a/r m and becomes negative for large a/r m , see e.g. the r m /£ = 1.0 
case in the left plot of Fig. [3j The phase transition, where the thermodynamic I = eigenvalue is a 



zero-mode, occurs precisely at the critical rotation identified in (2.7), namely at a/V m ^_ ~ 0.337 
(for r m /£ = 1.0 and in d = 6). This is what we expect from the thermodynamics discussion of 
subsection |2,2[ in regions B,C of Fig. [2] the thermodynamic I = eigenvalue is negative while in 
region A it is positive. The agreement between the numerical and thermodynamic results is an 
important check of the numerical code. 

The next interesting property of the spectrum of negative modes is encoded in the second curve 
counting from the left on the right plot of Fig. [3j According to Fig. [2j black holes with r m /£ < 0.864 
in d = 6 have a negative (positive) thermodynamic I = 1 eigenvalue if the rotation parameter is 



such that the black hole is in region C (region B), with the / = 1 zero-mode curve defined by (2.8). 
The black hole family displayed on the right plot of Fig. [3] obeys these conditions, since r m /£ = 0.5. 
Indeed, for rotations below (above) the critical rotation a/r m \ l _ 1 ~ 1.10, the thermodynamic I = 1 



eigenvalue is positive (negative) as predicted by (2.8). 
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Figure 3: Dimensionless negative modes Ar^ of the singly-spinning MP-AdS black hole in d = 6 
as a function of the dimensionless rotation parameter a/r m for fixed r m /L In the left plot, we 
represent the spectrum for r m jl = 1.0, while the right plot describes the spectrum for r m jt = 0.5. 
The latter has several (supposedly infinite) branches of zero-modes (for which A = 0) and the 
corresponding negative eigenvalues are labeled by the integer I. For this particular black hole 
family with r m /£ = 0.5, the values of a/r m at which the first few branches intersect the A = axis 
are: a/r m \ [=1 ~ 1.10, a/r m \ l=2 ~ 1.49, a/r m \ l=3 ~ 1.69, a/r m \ l=4 - 1-81, and a/r m \ l=5 ~ 1.88. 
As opposed to the thermodynamic Z = 0, 1 zero-modes (curve on left plot, first two curves on right 
plot), the zero-modes with I > 2 describe the onset of ultraspinning instabilities of the black hole. 



The I = and / = 1 negative eigenvalue curves are the only ones whose existence can be 



predicted by the analytical thermodynamic analysis of subsection 2.2 As discussed there, Refs. [381 



[2"7t [2"U] proved that a zero-mode of the thermodynamic Hessian (2.6) is also a zero-mode of the 



(modified) Lichnerowicz operator (3.5), but the latter can have additional zero- modes, unrelated to 
that Hessian. The thermodynamic instabilities associated with the / = 0, 1 negative modes do not 
correspond to classical gravitational instabilities of the black hole. In particular, the thermodynamic 
I = 0, 1 zero-modes do not describe the onset of a classical instability. Instead, they describe 
gravitational perturbations that change the mass and the angular momentum of the black hole, 
within the MP-AdS family, but preserve the temperature and the angular velocity, due to the 
boundary conditions. 

The thermodynamic / = 1 zero-mode curve defines the ultraspinning surface of the system. It 
was conjectured in Ref. [18] that only black holes which are outside this surface (i.e. in region C of 
Fig. [2]) might develop new zero/negative modes that have no thermodynamic interpretation, and 
that describe classical ultraspinning instabilities which are axisymmetric. The right plot of Fig. [3] 
confirms that the spectrum of MP-AdS black holes indeed has (a supposedly infinite sequence of) 
new branches of non-thermodynamic zeros- modes, for rotations strictly higher than the critical 



rotation a/r r< 



U=l 



associated with the I = 1 zero-mode. These families exhibit an underlying 
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harmonic structure (although the equations that we solve do not seem to separate into radial 
and angular equations). We use this property to suggestively label the several new branches by 
successive integers I = 2, 3, 4, . . .. This is an appropriate notation since the integer I coincides with 
the number of nodes that the metric perturbations have on the horizon y = |^] The values of 
o-l fm at which the first few branches intersect A = for the particular family of black holes with 
r m/f- = 0.5 in d = 6 are summarised in the caption of Fig. |3j The spectrum of black hole families 
with r m /£ < 0.906 is qualitatively similar to this particular case. 

In Fig. [4j we plot the number of negative modes detected numerically, for a grid of points in the 
parameter plane (r m /£,a/r m ) in d = 6 (left plot) and in d = 7 (right plot). This figure completes 
the information displayed in Fig. [2] (in particular, for d = 6, the black curves are the same plotted 
in Fig. [2]). Not only it includes the numerical and analytical results for the thermodynamic I = 0, 1 
modes but, in addition, it adds the numerical results relative to the I > 2 modes that have no 
thermodynamic interpretation. In region A (see Fig. g), we find no negative modes. In region 
B, the system has only the I = negative mode (blue dotted region in Fig. [4]). In region C, the 
/ = 1 negative mode is also present, but additional negative modes appear for higher rotations. 
In Fig. [4j these are signaled by the red, purple, green, yellow and pink dotted areas. In the red 
area, only the I = 0, 1 modes are negative, and there is no classical instability. However, for higher 
rotation, in the purple area and above it, the I = 2 harmonic becomes another negative mode of 
the modified Lichnerowicz operator (3.5). The transition from the red into the purple area defines 
the non-thermodynamic I = 2 zero-mode curve. This surface signals the onset of the ultraspinning 
instability in the parameter space. The boundary of the purple and green areas identifies the non- 
thermodynamic I = 3 zero-mode curve: above it the I = 3 harmonic is also excited, i.e. it becomes a 
negative mode. Similarly, the green/yellow (yellow/pink) transition marks the non-thermodynamic 
I = 4 (1 = 5) zero-mode curve above which the I = 4 (I = 5) harmonic is also a negative mode. Not 
shown in this figure, there should be an infinite sequence of new excited harmonics. These would 
become visible if we considered higher values of a/r m in region C. 

In Fig. [5j we proceed to a thorough analysis of the singly-spinning MP-AdS parameter space. 
We represent the data in Fig. [4] for d = 6 in terms of the dimensionless asymptotic conserved 
charges, M£ d ~ 3 and J£ d ~ 2 , instead of using r m and a to specify the black hole. We do this for the 
sake of clarity, because the charges are more physical, and also because the map (r m , a) — > (M, J) 
in Eqs. (2.2) is singular as \a\ — > £, so that the parameter space looks very different. The physical 
parameter range is strictly below the curve \J\ = Mi. To compare this figure with Fig. [ij-left, 
notice that: (i) the asymptotically flat limit, which was the r m /£ = line in Fig. Q-left, is the point 
(0, 0) in Fig. [5j (ii) the line a = £, for finite r m /£, in Fig. [Ij-left is now the asymptotic corner given 
by the limit R — > 00 of (R, R) in Fig. [5j so that the various coloured regions do not all meet in 
the parameter space as suggested in Fig. [4[ left; and (iii) the limit a/r m — > 00, within a < £, in 
Fig. [4]4eft is the curve J = M£ that puts a bound on the parameter space in Fig. [5] Notice that, for 
relatively small fixed values of the mass (in AdS units), say M/£ d ~ 3 = 1, there is always an I = 
negative mode, as in the asymptotically flat limit. On the other hand, there is a I = zero-mode 
for larger masses, say M/£ d ~ 3 = 2. For any fixed value of the mass, we find that a sufficiently high 



8 For the particular case where £ — > 00, this statement is illustrated in Fig. 3 of [19| . where we plot the functions 
5\(x, y) and 5cj(x, y) for I = 1 and I = 2. The figures are qualitatively similar for finite £. 
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d= 6 



d= 7 




Figure 4: Number of negative modes of the singly-spinning MP-AdS black hole in d = 6 (left) and 
d = 7 (right). The plots describe the dimensionless rotation parameter a/r m as a function of the 
dimensionless mass-radius r m /£. As we move from the bottom to the top, the new colored/dotted 
areas represent regions where a new negative mode of the (modified) Lichnerowicz operator gets 
excited: blue (one negative mode), red (2 n.m.), purple (3 n.m.), green (4 n.m.), yellow (5 n.m.), 
pink (6 n.m.), . . . The interpretation of the several curves plotted is described in Fig. [2] (Note that 
the faults in these figures, e.g. where we expected to find blue dots, correspond to parameter space 
points where our numerical code failed. They do not correspond to black holes with no negative 
modes.) 



rotation leads to the ultraspinning instabilities, in correspondence with Fig. [T] 

Our study is limited to stationary perturbations. Therefore, we have not proven that our 
I > 2 negative modes do indeed correspond to a region in parameter space where some black hole 
perturbations grow exponentially with time. Including time-dependence is not conceptually difficult 
but it is only computationally harder. Additional metric components would have to be perturbed 
(e.g. hf r , . . . ) in (3.1), and we would have to solve numerically a system of many more PDEs. 
This problem is however of fundamental interest since its solution would provide a definite proof of 
the ultraspinning instability in MP-AdS black holes, together with information on the instability 
timescale. Having said this, we nevertheless claim that we have found the stationary zero-modes 
that signal the onset of the ultraspinning instability. Our confidence comes from the fact that the 
analogous stability problem, including time-dependence, was studied in [20] for asymptotically flat 
cohomogeneity-1 MP black holes (equal angular momenta, odd d). There, the analogues of the 
I = 1, 2, 3, 4, . . . zero-modes are also present, and the corresponding perturbation sectors decouple. 
The time-dependent analysis confirmed the absence of a classical black hole instability in the I = 1 
sector of perturbations, and the existence of an ultraspinning instability in the I > 2 sectors. 
We take this result as good evidence in support of a similar interpretation in the singly-spinning 
MP-AdS case. 
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Figure 5: Number of negative modes of the singly-spinning MP-AdS black hole in d = 6. This 
figure has the data of Fig. now plotting the dimensionless angular momentum Jjl as a 

function of the dimensionless mass M/e d ~ 3 . Regular MP-AdS black holes exist for \J\ < ML We 
represent the location of the I = and the I = 1 zero-modes, for reference. 



The ultraspinning I > 2 zero-modes do not admit a thermodynamic interpretation, as we have 
stressed. They should correspond not only to the onset of classical instabilities, but also to the 
bifurcation to new stationary AdS black hole families with 'pinched' horizons, which have the same 
isometries as the MP-AdS black holes. This is the conjecture of Ref. |23j (following |14l HO] in the 
asymptotically flat case), represented in Fig. [TJ In particular, the I = 2 zero-mode should connect 
to a family of AdS black holes which are 'pinched' on the poles. The location of the I = 2 zero-mode 
is represented as point I in Fig. [TJ At some point along the phase diagram of this new family, a 
horizon topology phase transition should connect it to the AdS black ring with the same isometries, 
which has been constructed perturbatively in the limit of large dimensionless angular momentum 
|23j . Similarly, the I = 3 (Z = 4) zero- mode, identified as point II (III) in Fig. [TJ should connect to 
a new family that interpolates between the MP-AdS family and the AdS black Saturn (concentric 
rings), constructed perturbatively in [23 1 . Our results support the conjecture of Ref. [23] not only 
because we find the bifurcation points, but also because the shape of the numerical perturbations, 
with its harmonic structure labelled by I, is consistent with the appropriate 'pinches' of the horizon 
(footnote |J). 



As discussed in subsection 2.3 all the ultraspinning unstable MP-AdS black holes in the purple, 
green, yellow, pink, . . . dotted areas of Fig. [4] have Uh£ > 1, and thus all of them are also afflicted 
by the superradiant instability [29j |30j [31] . The instability is expected to be inherited by the 
'pinched' black hole families, at least close to the bifurcation point. 

An interesting open problem is what is the interpretation of the ultraspinning instability and 
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of the associated new AdS pinched black hole phases in the holographic dual field theory. Steps in 
this direction have already been taken in the context of the Scherk-Schwarz compactified AdS. In 
this case, in the hydrodynamic limit of the holographic theory one finds: i) pinched plasma balls 
|47] that should be in correspondence with the gravitational ultraspinning system, ii) deformed 
plasma tubes [38] that should be dual to the Gregory-Laflamme instability of Scherk-Schwarz- AdS 
black strings and iii) rotating plasma ball instabilities [39] that should describe the gravitational 
bar-mode instability \32\ 133] mentioned in the Introduction. 
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